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1. Introduction 



In the lattice QCD, most properties of hadrons are extracted from the hadronic correlation 
functions. The spectral function (SPF) has particular importance, since it may contain information 
beyond the stable ground state and a few excited states which can be extracted by standard fitting 
techniques. Recent development of analysis techniques such as the maximum entropy method 
(MEM) [|l|] have enabled direct extraction of SPFs from numerical data of lattice QCD simulation. 
At zero temperature, MEM has been successfully reproduced correct features of the SPFs [JlJ, ^. 

At finite temperature, we can calculate the SPF from the thermal green functions in principle 
using the same procedure as at zero temperature [|3|, ^. In particular, charmonium states have drawn 
much attention, since they probe the QCD plasma state through the changes of their properties 
[^, and hence are potential signal of the formation of quark gluon plasma in the heavy ion 
collision experiments |^. Several groups have studied the SPF of charmonium in finite temperature 
lattice QCD using MEM [||, ^ [T^ and their results indicate persistent //i/A state even above Tc. 
MEM has also been extensively applied to various areas of lattice field theories [jlT]]. 

While MEM is a powerful tool to extract SPF, it has intrinsic subtlety when applied to lattice 
QCD data of correlators. In this paper, we point out how each ingredient of MEM analysis causes 
such subtlety, focusing on an application to the correlators at finite temperature. In the next sec- 
tion, we consider general problems of MEM, and then in Sect. ^ describe particular problem at 
finite temperature caused by short extent in the temporal direction. Details of these analysis were 
presented in Ref. [||]. 



2. Maximum entropy method 
2.1 Outline of MEM 

First we briefly summarize the outline of MEM basically following Ref. [|l|], which reviews in 
detail MEM applied to data of lattice QCD simulation. We obtain the SPF, A{(o), from the given 
lattice result for the correlator, C{t), by solving the inverse problem, 

poo 

C{t)= dcoK{t,(o)A{(o), (2.1) 
Jo 

where the (continuum type) kernel K{t, ft)) is given by 

To extract the SPFA(ftj), MEM maximizes a functional Q{A;a) = aS[A] —L[A]. L[A] is the usual 
likelihood function, and minimized in the standard fit. The Shannon-Jaynes entropy S[A] is 
defined as 

'A{(oy 



poo 

S[A] = / dco 
Jo 



A{(o)-m{co) -A(ftj)log 



mi(0] 



(2.3) 



The function m{co) is called the default model function, and should be given as a plausible form of 
A{co). At the last stage of calculation the parameter a can be integrated out by a weighted average 
of prior probability for a. 
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Figure 1: Left panel: Samples of eigenfunction for the kernel, Eq. (2.2). The eigenfunction exp(M,) corre- 
sponds to a part of Eq. (2.5). Center and right panels: The results of mock data analysis. The dotted line is 
an original SPF and the solid line is reconstructed result by MEM. 



2.2 Singular Value Decomposition 

In the maximization step of Q{A;a) the singular value decomposition of the kernel K{t,(0) 
is usually used ^} Then the SPF is represented as a linear combination of the eigenf unctions of 
K{t,(o): 

A{(o) =m(»)exp< 'Y^biUi{(o) > , (2.4) 

where Ns is the number of eigenfunctions, bi are parameters, and m,(g)) the eigenfunction of the 
kernel K{t, (o). The number of degrees of freedom of A{(o) is accordingly reduced to the number 
of data points of the correlator. Although hi can in principle be determined uniquely from the data 
without introducing an entropy term, the small eigenvalues ofK{t,(o) lead to a singular behavior of 
the SPF; hence truncation of the terms is practically required, i.e. A^, may be less than the number 
of data points [[l3|]. In MEM, the entropy term stabilizes the problem and guarantees an unique 
solution for the coefficients of the eigenfunctions [|l|]. 



An outstanding feature of Eq.([2.4|) is that it can be fitted to generic shape without restriction 



to specific forms such as a sum of poles. However, the resolution of course depends on the number 



of degree of freedom in Eq. ([2.4|), and also on (O. An example of eigenfunctions is displayed in the 
left panel of Fig. |]. This figure indicates that the resolution of the function becomes worse in large 
(O region, because the superposed functions do not have enough variation. 

This feature of the eigenfunction is also shown in the mock data analysis. Center panel of 
Figure |I] shows the original (input) SPF and the reconstructed SPF from the correlator which is 
constructed by the original SPF with a random Gaussian noise. The original SPF has three peaks of 
the same width and hight at each co. When the noise of correlator is not so small, the reconstructed 
SPF does not agree with the original one; there is a tendency that the peak becomes broader than 
the original one at high energy region. The peak positions are correctly reproduced in this case. 

The right panel of Figure [I] shows more interesting example. When the width of a peak is 
narrow, MEM reproduces the shape rather well. However, for a case of large width, MEM fails to 
reproduce the shape of the original peak. 



Analysis of MEM without singular value decomposition was examined in Ref. [13]. 
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2.3 Default model function 



As mentioned in section 2. 1 , MEM needs a default model function to define the entropy term 
Q. Since small difference between a trial SPF and default model function makes the entropy term 
large, the default model function strongly affects the result of MEM when the quality of data is not 
sufficient. Therefore the default model function should include only reliable information we know 
beforehand. If not so, there is a risk the result might be controlled by hand. 

In the case of QCD, prior knowledge for SPFs is not so many, e.g. positivity and perturbative 
behavior at high energy region. In the case of point correlators, a natural choice of the default model 
function is the asymptotic behavior of the meson correlators at large CO in perturbation theory. We 
should remember, however, that such an asymptotic behavior is not observed in practical simulation 
because of the finite lattice cutoff. 

The risk caused by lack of reliable default model function is reduced by the quality of data. 
In fact, this perturbative form has been successfully applied to problems at zero temperature [jlj, |2|]. 
When we do not have good quality of data and reliable default model function, we have to check, 
at least, a default model function dependence to estimate an systematic uncertainty for the results. 



3. Application to finite temperature lattice QCD 

Since a temporal lattice extent, Nt, is restricted to l/Tut on finite temperature lattices, it is 
usually difficult to keep good quality of data as compared with zero temperature. Therefore we 
have to check the reliability of the results. 

One of good checks is to apply the method to the zero temperature data in the same condition 
as at finite temperature. To extract the SPF at finite temperature, at least, we should successfully 
reproduce the zero temperature SPF from the zero temperature correlator but the number of data 
points restricted to l/rt?,. We show these checks with our lattice data, which was obtained on an 
anisotropic lattice with j3 =6.10 and the renormalized anisotropy <^ = 4 having the spatial cutoff 
a^' =2.030(13) GeV[|]. 

We apply MEM to the correlator with restricted numbers of degree of freedom. The results 
with two types of such restrictions are displayed in Fig. ^ The left panel shows the dependence 
of the result on tmaxi the maximum t of the correlator used in the analysis. This case corresponds 
to the situation at T > 0. MEM fails to reproduce even the lowest peak for tmax < 16. The center 
panel shows the results when one alternatively skips several time slices in the analysis. This case 
corresponds to the coarsening of the temporal lattice spacing. Even for tg^p = 8 for which the 
number of data point is 6, MEM at least reproduces the correct lowest peak position while the 
resolution is not enough. These result indicate that the physical region of the correlator as well 
as the number of the degrees of freedom is important for MEM to work correctly. The required 
region of C{t) in the above analysis is tmax > O{0.5fm), which is not fulfilled around T ~ Tc. This 
situation may be improved by smeared operators. The left panel shows the results of MEM for the 
correlator with smeared operator. It is stable under the above two kinds of restriction for t,„ax of 
interest; at least the lowest peak position is correctly reproduced. 

Next we show the default model function dependence for a correlator with the smeared op- 
erator at finite temperature in Fig. ^ In this analysis we adopt the default model function of 
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Figure 2: Tests of MEM for restricting number of degree of freedom in the zero temperature correlator The 
first two panels show the results with a point operator, and the last is with a smeared operator. 
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Figure 3: The default model function dependence of SPF for a correlator with the smeared operator at finite 
temperature. 



m{(o) = mdm^o^, where mdm is determined by the perturbative asymptotic behavior of large ft) 
region. As mentioned in Sect. 2.3, since the default model function is not justified for lattice QCD 
simulation, we observe the default model function dependence in order to estimate a systematic un- 
certainty of the MEM results. In Figure ^ we change the mdm by factor 10, 1 and 0.1 respectively. 
The peak position is stable while the width is rather sensitive to such change of the parameter moM- 
These results indicate that it is difficult to discuss width of SPFs quantitatively based on MEM 
analysis. 



4. Conclusion 



In our previous paper ^ we have concluded the MEM is not sufficient for quantitative study 
of the SPF from our lattice data even if some smeared operators are adopted. If we roughly know 
the shape of SPFs, standard fit (or constrained curve fitting [14]) is rather appropriate for quanti- 
tative studies. Therefore we used MEM to find a rough image (fit-form) of SPFs and performed 
fit (or constrained curve fitting) with this functional form, such as multi Breit-Wignar type function, 
for more quantitative estimate of the width of the peak. In conclusion, although MEM is powerful 
tool to extract the SPFs from correlators, we have to use it carefully taking its subtlety into account. 
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The simulation has been done on NEC SX-5 at Research Center for Nuclear Physics, Osaka 
University and Hitachi SR8000 at KEK (High Energy Accelerator Research Organization). T.U 
has been supported by the JSPS Research Fellowships for Young Scientists. 
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